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1. Introduction 


The long-standing mystery in astronomy concerning the origin of the diffuse interstellar 
absorption bands (DIBs) was recently conveyed to Physics Today readers in an article 
{Huffmann, 1991) primarily concerned with describing the new field of fullerene research. 
Figure 1 (adapted from the above article) depicts the measured spectral dependence of the 
extinction that is generally attributed to interstellar dust. In Fig. la the optical density over a 
large energy range is shown. A smoothly varying main component is seen generally to in- 
crease with photon energy - displaying, however, an enormous "bump" at X ~ 2175A (5.7 eV). 
The diffuse interstellar bands appear as fine structure in Fig. la between about 1.9 and 2.8 eV. 
The DIB intensities generally correlate with the smoothly varying component of extinction 
("reddening") seen in a given direction. In Fig. 1b, the appearance of DIBs with the smoothly 
varying extinction component subtracted is depicted. 

A convincing explanation of the origin of the DIBs has evidently thus far not materialized. 
They are generally believed to be of interstellar origin - hence, the term: diffuse interstellar 
bands. Most recent explanations are solid-state based. Many postulate the existence in 
interstellar space of large organic molecules, such as polycyclic aromatic hydrocarbons 
(PAHs), long chain unsaturated hydrocarbons, fullerenes, etc., perhaps adsorbed onto various 
(unspecified) dust particles. With such models, serious difficulty is encountered in trying to 
account for the seemingly universal presence in interstellar space of the large and complex 
molecular structures whose existence is postulated. Furthermore, spectra predicted by the 
various models proposed thus far do not come very close to matching the fascinating and in- 
tricate spectral signature presented by the DIBs. 

in this paper, a new model is proposed to account for the DIBs. In this model, the DIBs 
result from a non-linear effect: resonantly-enhanced two-photon absorption of H£ ions located 
near the surface of the Stromgren sphere that surrounds an O- or B-type star. The strong light 
(at co i) that is required to "drive" the two-photon transition is provided by L(a) light emerging 
from the Stromgren sphere that bounds the H II region surrounding the star. In Sec. 2, a value 
—200 /zW/cm 2 is estimated for the L(a) flux at the Stromgren radius, R s , of a strong (05) star. 
(A pure hydrogen star is assumed.) It is shown in Sec. 3 that a c.w. L(a) flux of this intensity 
should be sufficient to induce a few percent absorption for visible light radiated by the same 
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star at a frequency {a> z ) that completes an allowed two-photon transition, provided ( 1 ) the 
L(«) radiation happens to be nearly resonant with the frequency of a fully-allowed absorber 
transition that effectively represents the first step in the two-photon transition, and ( 2 ) an ef- 
fective column density ~10 18 /cm 2 of the absorber is present near the Stromgren sphere radius, 
R s . In Sec. 4, it is shown that a few almost exact resonances (i.e., with frequency offsets from 
L(a), A, not exceeding ~2 cm' 1 ) are found in the allowed Hi system: 2 p 7 r l/ (v', J') 4 - 1s<r g (v", 
J"). It is then indicated that the "second step" Hi transitions: [4d cr g , 5ga^, 3dr5 5 , 5gn g , Adn g , and 
3s«r ? ] 4 - 2pn u (v\ J') can possibly account for the DIBs. This visible absorption spectrum is, 
in principle, exactly calculable. To begin to explore whether or not there might really exist 
exact coincidences between measured locations of DIBs and transition wavelengths predicted 
by the model, one could initially attempt to compute vib-rot level energies for the (four) above 
even-parity states that are bound, preferably with adiabatic corrections incorporated into the 
Born-Oppenheimer approximation. 

The required presence of Hi molecules located near R s could conceivably result from di- 
rect photoionization of H 2 molecules located in a compact shell that serves as the transition 
region between the H II sphere and the cold surrounding hydrogen gas. In this case, one 
would expect the Hi ions to be produced in the ground state with a vibrational distribution 
given by the known Franck-Condon factors for ionization of hydrogen. It is, therefore, quite 
interesting that the linear Hi absorption spectrum one calculates for such a distribution has 
generally similar features to those possessed by the smoothly-varying component of the 
measured "interstellar" extinction curve, including the large "bump" at ~2175 A. 


2. L(a) radiation emerging from radiativeiy excited nebulae of pure hydrogen (H II regions) 

The source of excitation of the gas is assumed to be an embedded bright star of early 
spectral type - for example, an O- or a B-type star. We rely here on the simplest model for this 
gas - as described, for example, in Chapter 5 of the text by Dyson and Williams (Dyson and 
Williams, 1980). Table 1, adapted from Table 5.3 of one of Spitzer's books (Spitzer, 1978), 
shows what the total rate of photon emission shortward of the Lyman limit for typical O and 
B stars is. We here assume that our embedded star is like the 05 one listed, with N u ~ 5 x 
10 49 s- 1 . Note that the radius of the star is only 13.8 x R 0 , or 3.2 x 10 " 7 parsec (pc). 


Spectral Type 

T C ff(°K) 

R/R q 


05 

47,000 

13.8 

51 x 10 48 

06 

42,000 

11.5 

17.4 

07 

38,500 

9.6 

7.2 

08 

36,500 

8.5 

3.9 

09 

34,500 

7.9 

2.1 

BO 

30,900 

7.6 

0.43 

B1 

22,600 

6.2 

0.0033 


Table 1. Radiative Properties of Early-Sequence Stars 
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Let us estimate the L{«) flux emerging from the "Stromgren sphere" (the sphere whose 
radius, R s , delineates the boundary between the almost fully ionized H II interior, and the cold 
hydrogen gas exterior). We have, approximately, 

47r/?!'F L ( a) = y Tt R l n lp 2 {T e ) (1) 

This equation assumes that there is one L(a) photon produced for every recombination (to 
n>2) that occurs. Also, if a L(a) photon is absorbed by an H atom, it will simply be re-scattered 
in another direction. All L(«) photons produced will eventually escape through the surface of 
the Stromgren sphere. 


An excellent numerical approximation for the total recombination coefficient, /? 2 (T C ), is 
(Dyson and Williams, 1980): 

^(T e ) = 2x10- 16 r- 3/4 m V 1 , (2) 

where T c is the temperature of the gas in the H I! region. Estimates of such temperatures lie 
roughly in the range 5,000 - 10,000 K, with rather little variation from one nebula to the other 
(Dyson and Williams, 1980). We assume T c = 10,000 K for our case. Then /? 2 ~ 2 x 10" 19 m 3 
s~\ 


The Stromgren radius, R s , is given by (Dyson and Williams, 1980): 


3 N u 1/3 

r s = (—T~) m. 


47TD e /?2 


Thus, combining this with Eq(1), one has: 

L(«) = ~3 


y - 1 n 4/3 /? 2/3 n N(J ) 3 

T L(«) ~ v n e P 2 W An ) • 


(3) 


(4) 


Let's assume that n c - 10 5 cm- 3 . Then one has: = 1.21 x 10 18 photon s- 1 m~ 2 . This is a 

L(a) intensity of about 200 /xW/cm 2 . 


For the parameters chosen, the value of R s is about 1.8 x 10 15 m, or about 0.06 pc. - some 
2 x 10 s times greater than the radius of the embedded 05 star. If one multiplies 47 t times the 
square of R s times the L(a) flux calculated above, one gets back, of course, 5 x 10 49 photons 
s _1 . 


From the discussion given above, the linewidth of the L(a) radiation emerging from an H 
II region has to be the Doppler width at the temperature T e . This is given by the formulas 
(Spitzer, 1978): 


with 


Av = 


2( In 2 ) 1/2 b 
X 


y ^ /2 

b = 1.290x1 0 4 (-~) cm/ sec. 


(5) 

( 6 ) 
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Choosing T e to be 10,000K, A (atomic weight) to be 1, and X to be 0.1215 x 10~ 4 cm gives Av 
cm- 1 . Observed DIB linewidths can be both larger and smaller than the 1(a) linewidth 
assumed in our model. 


3. Resonant enhancement of two-photon absorption cross-section 

For linearly polarized light waves, the absorption cross-section of an atom (or molecule) 
for light at frequency v 2 induced by light at v 1 for transitions from the ground state g to the 
excited state f is given by (Bjorkholm and Liao, 1974): 

2 

a(v 2 ) = 1. 15x10“ 34 |^T < f\z\n >< n\z\g > [(£„ - hv 2 )~ 1 + ( E n ~ hv \Y ’ll v 2p( v 2 + v iKi crr)2 ' ( 7 ) 

n 

where /•, is the intensity in W/cm 2 of the radiation field at v-,, p(v 2 + v,) is the normalized 
lineshape function for the (Hi) transition, the matrix elements are expressed in Bohr radii, and 
the summation is carried out over all intermediate states of energy E n , expressed in rydbergs. 
For the 3s-4d transitions in Na (with the 3p levels as the dominant intermediate states), use 
of Eq. (7) leads to the result (Bjorkholm and Liao, 1974) that the induced cross-section for this 
transition is approximately (5 x 10" 16 ) h cm 2 when E 3p - hv., = 1.0 cm" 1 . With an effective de- 
tuning of 1.0 cm -1 and with /-, = 200 /A/V/cm 2 , a column of (for example, HI) ground state 
molecules of length 10 16 cm (0.003 pc) and of density 100/cm 3 could thus, in principle, absorb 
10% of the visible background light at those frequencies completing two-photon resonances. 
Note that the total HI column length (0.003 pc) needed here is only about 5% of the Stromgren 
radius (0.06 pc) in this case. 


4. HI transitions that are nearly resonant with L(a) and those that complete the two-photon 
absorption process 

If a vibrational level v" of the 1sa g ground electronic state has a good Franck-Condon 
overlap with a vibrational level v' of the 2pn u excited state, transitions between the two levels 
are strongly allowed. Using adiabatically corrected values (Bishop et al., 1975) for the vib-rot 
levels of these two electronic states, one can show that the P(1) line of v' = 5 <- v" = 15 almost 
exactly coincides with the peak wavelength of the L(a) radiation, and that the P(1), Q(2), and 
R(3) lines of v' = 11 <- v" = 18 are all within 2 cm -1 of the L(a) peak. The DIB wavelengths in 
the model proposed here should thus correspond to wavelengths of allowed transitions from 
(low J' value ) 2pn u states with v'=5 and v'—18 to even-parity states dissociating to H + + 
H(n>3). Although Born-Oppenheimer potentials for the twelve states dissociating to H + + 
H(n = 3) have been calculated (Madsen and Peek, 1971), vib-rot level energies for the (bound) 
states have not, for example, been determined. Hence, until accurate calculations for these 
states become available, the best procedure for estimating DIB wavelengths from the model 
is to determine vertical energy separations between upper and lower potential curves at the 
right-hand turning points of v' =5 and v' = //, i.e., at 13.0 a 0 and at 20.5 a 0 (Fig. 3). We illustrate 
by carrying out this procedure for the 3dr5 ff <- 2pn u transition. This has a large, constant (i.e., 
as a function of R) dipole moment of 3.00 a.u. (Ramaker and Peek, 1972). A vertical line drawn 
at 20.5 a 0 intersects the 3d<5 g potential curve at a point ~282 cm -1 below its dissociation en- 
ergy. The upper level for this transition is thus stable and would not contribute to the DIB 
linewidth, which would in this case presumably be primarily determined by the spectral extent 
over which resonance enhancement given by Eq. (7) occurs. The predicted DIB wavelength is 
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here 6679 A. Could this possibly be the strong, narrow (width at half depth ~ 2.5 cm -1 ) DIB 
observed at 6614 A? At 13.0 a 0 , a vertical line intersects the 3d<5 ? potential at a point ~178 
cm -1 above its dissociation energy (and 125 cm -1 above the height of the shallow "barrier" 
existing at R = 37.5 a 0 ). The DIB transition, predicted in this case to occur at 6156 A, should 
thus be strongly broadened. Could this be the strongly broadened DIB (width at half depth ~ 
80 cm" 1 ) that is observed at 6177 A? In Table 2, we show the wavelengths that are generated 
in carrying out this procedure for all six even-parity states dissociating to H + + H(n = 3). 
Wavelengths predicted to be sharp (s) or diffuse (d) are so indicated. 


Upper State 

R - 13.0 a 0 

R = 20.5 a 0 

5 9*sr 

4054 A (d) 

5959 A (d) 

4d n g 

4658 A (d) 

5749 A (d) 

366* 

6156 A (d) 

6679 A (s) 

5ga g 

4669 A (d) 

9219 A (s) 

4dc Tg 

6215 A (s) 

6842 A (s) 

3s <jg 

4467 A (d) 

5441 A (d) 


Table 2. Predicted DIB Wavelengths 


5. Is Hi possibly also the source of the "bump" at l n 2175 A? 

As stated in the Introduction, the required presence of Hi molecules located near R s could 
conceivably result from direct photoionization of H z molecules in a compact shell that serves 
as the transition region between the H II sphere and the cold surrounding (H 2 ) gas. In this case, 
one would expect the Hi ions to be produced in the ground state with a vibrational distribution 
given by the known Franck-Condon factors for ionization of hydrogen, ft is apparent from the 
Dunn and Oksjuk results (Fig. 4) that the cross-section for 2pov <- 1sa ff photodissociation for 
such a vibrational distribution peaks somewhere in the 2100 - 2200 A range. Thus, an intrigu- 
ing additional feature appears in the model presented here! 

One can even "push one's luck" and try to ask why the interstellar extinction curve re- 
verses its decline at ~8 eV, after it has fallen considerably from its 2175 A peak, and strongly 
starts to increase again with increasing photon energy. A possible explanation is that at 
shorter wavelengths the vibrational Hi 1s<7 3 distribution that results from H z photoionization 
starts to absorb again when the photon energy becomes high enough for the 2pn u 1sct ? 
transition to occur. There are also allowed transitions from 1s<7 ? to higher states of u parity, 
states dissociating both to n = 2 and higher n. 
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Fig. 1. Extinction due to interstellar dust (Huffman, 1991) 




Fig. 2. Some potential curves for Ht 
(Levine, 1974). 



Fig. 4. <t(2) [2pa u «- 1scr g transition] for the 
Ht vibrational distribution produced by 
photoionization of H z (Dunn, 1968). 


Fig. 3. 1sa ff and 2pn u potential curves (Beckel, C.L.; 
Shafi, M.; and Peek, J.M.: 1973, J. Chem. Phys. 59, 5288). 
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